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Abstract 

We demonstrate that electric dipole moments (EDMs) strongly constrain possi- 
^ , ble SUSY contributions to the CP asymmetries of B processes; LL and/or RR 

•r:^ I flavour mixings between second and third generations are severely restricted 

^ I by the experimental limit on the mercury EDM, and so therefore are their pos- 

CN ! sible contributions to the CP asymmetries of 5 — > (pK and B rj'K. We find 

^ I that SUSY models with dominant LR and RL mixing through non-universal 

O I A-terms is the only way to accommodate the apparent deviation of CP asym- 

[ metrics from those expected in the Standard Model without confiicting with 

I I the EDM bounds or with any other experimental results. 

(D : 

'f: : 1 Introduction 

> . 
•1—1 , 

^ ■ The most recent results of BaBar and Belle collaborations [1] on the mixing-induced 

' asymmetries of 5 — > (j)K and B rj'K indicate possible deviation from the Standard 
Model (SM) expectations. The Belle experimental values of these asymmetries are given 
by 

S^K = 0.06 ±0.33 ±0.09, (1) 

S^^K = 0.65 ±0.18 ±0.04. (2) 



The BaBar experimental values are 



S^K = 0.50 ± 0.25^°:°!, (3) 
S.n'K = 0.27 ±0.14 ±0.03. (4) 



Comparison with the world average CP asymmetry Sj/^k = 0.726 ± 0.03 shows that the 

'.11; 
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average CP asymmetry of all 6 — >^ s penguin modes from the Belle results is 0.43l^o;^^ 



which is 2.4(7 away from the SM result, and from the BaBar result is 0.42 ± 0.10, a 2.1a 
deviation. 

Supersymmetry (SUSY) is one of the most popular candidates for physics beyond the 
SM, and a natural place to look for explanations of such deviation. Indeed in SUSY 
models there are many new sources of CP violation besides the CKM phase. However 
stringent constraints on these phases are usually obtained from the experimental bounds 
on the electric dipole moment (EDM) of the neutron, electron and mercury atom. Because 
of this it is it a challenge for SUSY models to give a new source of CP violation that can 
explain the possible discrepancy between CP asymmetry measurements and the expected 
SM results, whilst at the same time avoiding the overproduction of EDMs. 

It is known [2-4] that SUSY models with a large squark mixing and order one phase 
between the second and third generations can accommodate the CP asymmetry results 
via gluino exchange. The squark mixings can be classified, according to the chiralities 
of their quark superpartners, into left-handed or right-handed (L or R) squark mixing. 
The left-handed mixings for the down-squark are given by the mass insertions {Sfi)ij 
and {Sfji)ij, and the right-handed mixings by {Sji]^)ij and {Sjiji)ij. It is remarkable that 
in order simultaneously to satisfy the measurements of S^px and S'^'k and explain the 
deviation between them, both left- and right-handed contributions have to be involved. 
This is because the left- and right-handed contributions have an opposite sign due to the 
different parity in the final states of B ^ (j)K and B rj'K [4]. 

In this paper we argue that a large flavour mixing between the second and third 
generation via (5^/,)23 and {6f^ji)23 leads to a large (5^^)22, which produces a large strange 
EDM and consequently overproduces neutron EDM (assuming the Parton model) and 
mercury EDM. We will show that, taking EDM constraints into account, the possible 
solution of the S^f^K and Sn'K discrepancy based on (5^2.)23 and (5^/j)23 is disfavoured. This 
leaves the scenario with large mass insertions {S'[ji)23 and ((^^i)23 (due to non-universal 
trilinear A-terms) as the only possible consistent solution. 

This paper is organized as follows. In section 2 we introduce the super symmetric 
contributions to the strange quark EDM which could be enhanced by large mixing between 
the second and third generation and leads to a large Hg EDM. Section 3 is devoted to 
imposing the EDM constraints on the SUSY phases in a model independent analysis, 
and the impact of these constraints on the SUSY contribution to the CP asymmetries of 
B — > (l)K and B — > rj'K. In section 4 we give numerical results and show correlations 
among the Hg EDM and the CP asymmetries of B-decays. Our conclusions are given in 
section 5. 

2 Supersymmetric contributions to strange quark EDM 

As mentioned in the introduction, SUSY models have several possible sources of CP 
violation in addition to the CKM phase. These CP phases can have important implications 
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for CP violating phenomenology. In particular they can induce large EDMs of quarks 
and leptons at the one-loop level that far exceed the experimental limits, and stringent 
constraints on SUSY CP phases are found [5]. The most recent measurements for the 
neutron {dn) and mercury (dng) EDMs lead to the following limits: 

dn = 6.3 X lO^^^ecm, (5) 
dn, = 2.1 X IQ-^^ecm . (6) 

The neutron EDM receives contributions of different sources and the predicted value in 
any particular model depends quite strongly on the particular model of the neutron used 
for the calculation. Because of this it is worth briefly summarizing the calculation. 

The major contributions to the EDMs come from electric and chromoelectric dipole 
operators and the Weinberg three-gluon operator: 

jC = -'-dfqa,,^,qF^^ - '-d^qa.^^.T^qG^^'' - ^d^/„,eG„^,G^,GeA.6^^^^ (7) 

In order to evaluate the neutron EDM, one needs to make some assumptions about the 
internal structure of the neutron. The models can be classified as follows; 

1- The chiral quark model 

In this model the neutron EDM is related to the EDMs of the valence quarks; 

4 1 

dn = gC?d - gO?«- (8) 

The quark EDMs are given by 

d, = ,''d^ + n''^d^ + ,"'^^, (9) 

where the QCD correction factors are given by rj^ — 1.53, t;*^ ~ r^*^ ~ 3.4 and where 
A ~ 1.19 GeV is the chiral symmetry breaking scale. 

2- The parton quark model 

Here one assumes that the quark contributions to neutron EDM are weighted by the 
factor defined as (w||?7^75?|w) = ^qSf^, where 5*^ is the neutron spin; 

dn = {^dd^ + A.rff + A,df ) , (10) 

where the individual quark contributions are given in terms of the gluino, chargiono and 
neutralino contributions 

dg = dl + dl^ +df. (11) 

3 



The following values for are usually used: A^^ = 0.746, A^^ = —0.508, and A^ = —0.226. 
The main difference between the parton quark model and the chiral quark model is the 
large strange quark contribution in parton model. Also in this model, the relevant con- 
tribution is only due to the electric EDM of the quarks in contrast with the chiral quark 
model where the chromoelectric and three-gluon operators contribute as well. 

3- QCD sum rules 

The QCD sum rules analysis of ref. [6] leads to the following relation between the 
neutron EDM and the electric EDMs and chromoelectric EDMs of u and d quarks: 

dn = 0.7(df - 0.25df ) + 0.55e(d^ + 0.5d^), (12) 

where the value of quark vacuum condensate {qq) = (225GeV)^ has been used. It can be 
seen from the above equation that the QCD sum rules cannot incorporate the effect of 
the strange quark in the neutron EDM. 

4- The Chiral Lagrangian approach 

In ref. [7], the chiral lagrangian approach was adopted to try to incorporate the strange 
quark chromoelectric EDM contribution to the neutron EDM. This analysis leads to the 
following result for the neutron EDM in terms of the quark chromoelectric EDM: 

4 = (1.6c?^ + 1.3(i^ + 0.26c?f)ecm. (13) 

Passing to the mercury atom EDM, the major contribution here comes from T-odd 
nuclear forces in 7r° and 7] couplings to the nucleus, which is generated by the chromoelec- 
tric EDMs of the constituent quarks. The resulting EDM of the mercury atom is given 
by ref. [8] as 

dH^ = -e(rf^ -d^- 0.012rff ) X 3.2 x lO'l (14) 

Although the coefficient for the d'^ is much smaller than the coefficients of the chromo- 
electric EDM of the down and up quarks, this contribution is still important since rf^ 
itself is enhanced by the heavy strange quark mass and by the relatively large mixing 
in the second generation. Recently the mercury EDM has been reconsidered in the light 
of the QCD sum rule calculations, with the result that the coefficients multiplying the 
first generation quarks could be reduced by a factor 2.5-3 [6] (see ref. [9]) for a recent 
discussion). Our study will depend mainly on the strange quark EDM so this uncertainty 
will not effect our conclusions. We will therefore use the older bound for this study, and 
comment at the end. 

The dominant 1-loop gluino contribution to the EDMs is given by 

= -l^Q^,u^lM&nM^{x), (15) 

O TV III ~ 

a 
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df = -l^Qs^lm{Si^),,M,{x), (16) 

6 n rri r 

a 

- 1^ 3 M5L)22 M,(x), (17) 

d 

where x — mymg. The current experimental bounds using the parton model for the 
neutron EDM imply the following constraints on the relevant mass insertions [5] : 

Im(5^^)n < 1.9 X 10-^ Im(5^^)22 < 6.6 x 10-^ (18) 

where to illustrate we have taken mj ~ 500 GeV and x — 1. The experimental limit on 
the mercury EDM leads to a stronger bound on the imaginary part of ((^iij)ii and about 
the same bound on the imaginary part of {S'lji)22- 

Im(5^^)n < 6.7 x 10-^ Im(5^K)22 < 5.6 x 10"^ (19) 

As alluded to above, the mass insertion (5^^)22 is more sensitive to the mixing between 
the second and the third generation, so the bound on its imaginary part is the relevant 
one for our analysis. It is remarkable that the bounds obtained on this quantity from the 
mercury EDM and neutron EDM are almost the same, however, as we emphasized, the 
computation of the neutron EDM is more model dependent. Therefore in our analysis we 
will concentrate on the constraint obtained from the mercury EDM. 

The explicit dependencies of {S'[^)22 and (5^^)22 on the LL and RR mixing between 
the second and the third generations are given by 



(<^L)22 = {51^)23 {Stn)33 {SU)32 + 5^23 (5^)33 (^^^32 , (20) 



f:d \ f xd \ I s:d 



where (5^^)33 = {S%l)33 ~ "'''^^''"f'"^^ Recall that the EDM is proportional to the 

d 

imaginary part of the coefficients of the d*i^dR term in the Lagrangian. In the MSSM the 
relevant part of the Lagrangian is given by £ ~ {YA — /xtan/3) d^dn + h.c. where h.c. 

refers to {YA - //tan/3)* dLd*R. 

The EDM imposes stringent constraints on the flavour conserving CP phases of the 
and /i terms. It is reasonable therefore to assume that these phases are suppressed, 
in which case (5^/^)33 ~ mb/mj ~ 10^^. Also, due to the hermiticity of the LL and RR 
sectors of the squark mass matrix, {5'Il{rr))^2 = {^ll{rr))23- Thus one finds 

{5i^)22 ^ 10-^[(<5^j23(O;3+((4«)23(<^fj;3)1 (21) 

Furthermore, {5'1^)22 can also be expressed as 

(^Lij)22 = (5^l)23 (^m)32 + [(4i?)23 (5rl)32] , (22) 

where (5^^)32 = (^^l)23- l^ ney± section, we will determine the values of lm.{5j^j^22 
and Im(5^^)22 within the regions of the parameter space that satisfy the experimental 
results of S^k and SrfK- We will show that the EDM of the strange quark allows the 
possibility of SUSY models with large LR (RL) mixing only. 



3 SUSY contributions to the CP asymmetries S(j,K 

and Sj^'K 

Including the SUSY contribution, the effective Hamiltonian H^^^ for these processes 
can be expressed via the Operator Product Expansion (OPE) as 



AB=1 
eff 



G p 
V2 



Ap I CiQi + C2Q2 + X] ^iQi + CllQ^l + CsgQsg 1 + H.c. 
p=u,c \ j=3 / 



(23) 



where \p = VpbVpg, with Vpb the unitary CKM matrix elements satisfying At + A„ + Ac = 0, 
and Ci = Ci{iJ,b) are the Wilson coefficients at low energy scale ~ nib. 

As emphasized in refs.[2,4], the dominant gluino contributions are due to the QCD 
penguin diagrams and chromo-magnetic dipole operators. At the first order in MIA, the 
the gluino contributions to the corresponding Wilson coefficients at the SUSY scale are 
given by 



Ci 



r<9 



2V2GFm1^ ^^^23 



V2G 



1 



--B^{x) + ^-B^ix) + ^-P,{x) + ^-P^ix) 
y^i(^) + ^B2{x) - ^P,{x) - ^P2{x) 
~B,{x) + ^-B2{x) + ^P,{x) + ^P2{x) 



m-r, /I 



{^II),^ [-^M,{x) + ^M,{x)) + {5i^)^^-^^ [-M,{x)+^M,{x 



23mb V3 



.(24) 



where Ci,8g are obtained from Cj^sg by exchanging L <-> i? in (5^^)23- It is clear that 
the part proportional to LR mass insertions in Cf^ which is enhanced by a factor mg/rrib 
would give a dominant contribution. Using the QCD factorization mechanism to evaluate 
the matrix elements, the decay amplitude of S — > (f)K can be presented in terms of the 
relevant Wilson coefficients as follows [2]: 



A{B ^ <I>K) = -t^mlF^-^ U E ^i(0)(Q + Q), 



(25) 



where Hi{(j)) arc given in Ref.[2] and the Wilson coefficients Cj and Cj are defined ac- 
cording to the parametrization of the effective Hamiltonian in Eq.(23) 



G 



(26) 



Therefore, the contributions of the RR and RL terms in R^f, have the same sign as the 
LL and LR ones. For instance, with rUq = rrig = 500 GeV, one obtains 



R, -0.14e-^°-^(<5iJ,3-127e-^°°«(5^^)23-0.14e-°-^(4^)23-127e-^°°«(4J,3. (27) 

From this result, it is clear that the largest SUSY effect is provided by the gluino con- 
tribution to the chromomagnetic operator which is proportional to (5ii^)23- However, the 
6 — > S7 constraints play a crucial role in this case. For the above SUSY configurations, the 
6 — > S7 decay constrains the possible gluino contributions since it sets |(5f^)23| < 0.016. 
Despite this, on implementing the bound in Eq.(27), we see that the gluino contribution 
(proportional to (^^7^)23) is still able to generate large values for R^, consequently driving 
S<pK towards the region of small values. 

Although B (f)K and B rj'K are very similar processes, the parity of the final 
states can vary the result. \n B (j)K the contributions from Ci and C-i to the decay 
amplitude are identically the same (with the same sign), whereas in i? — > if]'K they have 
opposite signs. This can be easily understood by noting that 

{(t>K\Qi\B) = {(t>K\Qi\B) . (28) 

which is due to the invariance of strong interactions under parity transformations, and 
to the fact that initial and final states have the same parity. However, in the case of the 
B — > r]'K transition, where the initial and final states have opposite parity, we have 

{'n'K\Q,\B)QCDF = -WK\Qi\B)QCDF. (29) 

As a result, the signs of the Cj and Ci in the decay amplitude are different for B — > rj'K, 
and so the sign of the RR and RL in i?^/ are different from the sign of LL and LR in 
contrast with the i?^ case. Using the same SUSY inputs adopted in Eq. (27), we have 

R,, -0me^'-'\Sir.)23 - 64(5^^)23 + 0.07e^°-2^(4«)23 + 64(4^23 (30) 

Following the parametrization of the SM and SUSY amplitudes in Ref . [4] , S^k can be 
written as 

^ _ sin 2(3 + 2R^(^n') cos 612 sin(g^(^.) + 2(3) + sin(2g^(^/) + 2(3) 

"^^^ 1 + 2-R^(^') cos S12 cos 9^(^r,') + Rl{r,') 

where R^ = |A^usy/^sm|^ ^ arg(A^usY/^sM)^ ^-^e strong phase. In order to 

accommodate the experimental results of S^k and Sr,'K we should have at least one of the 
following two scenarios [2,4]: large mixing between the second and the third generations 
in LL and RR sectors or large mixing between the second and the third generations in 

LR and RL sectors. 

As can be seen from Eq.(31), the deviation of S^(.ri')K from sin 2/5 strongly depends 
on the size of R^p^r]')- The minimum values of S^k and Sr^'x can be obtained by large 
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Figure 1: Hg EDM, S^k and Sn'K as function of the magnitude of (5^2.)23 (l^ft plot) and 
(5^^) 23 (right plot). The phases of the mass insertions are assumed to me 0{ti /2). Also 
nig = rriq = 400 GeV is used. 

values of \{S'[l)23\ ~ C(l), |('5^,k)23| ~ 0{1) and phases of (5^^)23 and {S'^r)23 of order 
one or 1(5^^)231 ~ 0(10-3), |(4^)23| ~ 0(10-3) and phases of (5^^)23 and (4J23 of 
order one. It is important to note that to have deviation between S'^x and Sn'K, the 
contributions from (5^^)23 (('^Li?)23) and (^hj?)23 (^{^rl)'2s) should be different, so that 
the gluino contribution to S^k becomes larger than its contribution to Srj'K [4]. It is also 
worth mentioning that, as can be seen from Eq.(31), the effects of LL and RR mixing on 
S(f){ri')K remain limited compared to the effect of LR and RL. 



4 EDM constraints on Sfpx cind S^^'x 

We now come to the main point of this letter, which is that such large values for the 
magnitudes and the phases of (^^2^)23 and (5r/j)23 may significantly enhance the strange 
quark EDM thereby overproducing mercury and possibly neutron EDMs. It is interesting 
to ask therefore whether SUSY is still able to accommodate such large magnitudes and 
phases, and if so, are they restricted. 

As mentioned in the previous section, there are two possible sources of enhancement: 
the first is the combined effect of (5^^)23 and (5^/j)23, the second source is either (5^j^)23 
or (5^/2)23 combining with (5^/j)32 or (5^^)32- However, within minimal flavour models 
such as minimal supergravity (where the trilinear couplings are universal), the size of the 
mass insertions (5^/2)23 and (5^^)23 are of order 10-^ and 10^'' respectively. Therefore 
the imaginary part of the induced mass insertion (5^^)22 can easily be below the bound 
obtained from the experimental limit on Hg-EDM. In Fig. 1, we plot both S^k, S^'k 
and dHg/{dHg)Exp as functions of \{5ii)22\ and \{5%j^)2i,\- We assume that arg[((5f^)23] ~ 
arg[((5^^)23] ~ 7r/2, in order to enhance their effects on the CP asymmetries of 5-decays. 
We also fixed mg = = 400 GeV. 

As can be seen from the figure, in these scenarios the values of Hg EDM are well below 
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Figure 2: Scf^K (left plot) and Sj^'k (right plot) versus the ratio of the Hg EDM to its 
experimental value. 

the current experimental limit. However, we cannot account for the CP asymmetries S^k 
and S-r^iK^ particularly Sr^'K which has been the subject of recent measurements by the 
BaBar collaboration. In this class of models with dominant (5i2,)23 or (5^/2)23 mass 
insertions, the value of S-^^'k is close to the SM prediction of sin 2/3. Therefore if the 
present S^'k result is confirmed, these models will be disfavoured. 

In Fig. 2, we display scattering plots for Stpx and Sr^iK versus the ratio of Hg EDM 
to its experimental hmit. We set mo — mg = 400 GeV. The other relevant parameters 
are scanned as follows: \{5'[i)23\ varies from to 1, arg[((5f^)23] and arg[((5^^)23] are in 
the region [— tt, tt]. As can be easily seen from Fig. 2, within the region of the parameter 
space where S'<^x and Sj^k fit the experimental data, the Hg EDM exceeds with many 
order of magnitudes its experimental bound. This imposes severe constraints on this 
scenario of simultaneous contribution from LL and RR mixing to accommodate both 
the experimental results of S^px and Sr^K- This result is in agreement with that of of 
Ref. [10]. Returning to the question of the precise numbers in the bound, it is clear from 
the figures that even if the strange quark contribution to the mercury EDM were reduced 
by a whole order of magnitude (rather than the factor 2.5-3 reduction implied for the first 
generation contributions to the Hg-EDM from the sum-rule calculations) this conclusion 
is unchanged. 

Therefore, we may safely conclude that SUSY models with dominant LL and/or RR 
large mixing between second and third generations will be ruled out if the experimental 
results of S^k, S^^k are confirmed. SUSY models with dominant LR and/or RL mixing 
via non-universal A-terms, seem to be the simplest way to account for CP asymmetry 
S^K and SriK without conflicting with EDMs or any other experimental results. 

Before we conclude, we give a quantitative prediction for the Hg-EDM due to the 
effect of large {5'lj^2z- As shown in Ref. [2], in order to accommodate the experimental 
result of the CP asymmetry |(5fij)23| should be of order lO""^ and arg[(5f^)23] ~ 7r/3. 
These values lead to S(j,K — 0.2. Assuming (the minimal assumption) that the soft scalar 
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masses are universal at the SUSY breaking scale, one finds that at the electroweak scale 
(^11)23 is of order 10^"^ and (5^/2)23 ~ 10~^. Hence one finds that (5^^)22 ^ 10^^ which 
implies that dng ~ 0-'2{dHg)Exp- We find it intriguing that the Hg EDM experiment is so 
close to testing CP violation in the flavour changing sector. 
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